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THE PROPER COMBINATION OF LIFT LOADINGS FOR LEAST DRAG 

ON A SUPERSONIC WING 1 


By Frederick C. Grant 


SUMMARY 

La, grange’s method of undetermined multipliers is applied to 
the problem of properly combining lift loadings for the least 
drag at a given lift on supersonic wings. The method shows 
the interference drag between the optimum loading and any 
loading at the same lift coefficient to be constant. This is an 
integral form of the criterion established by Robert T. Jones for 
optimum loadings. 

The best combination of four loadings on a delta wing with 
subsonic leading edges is calculated as a numerical example. 
The loadings considered have finite pressures everywhere on the 
plan form. Through the sweepback range the optimum com- 
bination of the four nonsingular loadings has about the same 
drag coeffi.cien t as a fiat plate with leading-edge thrust. 

INTRODUCTION 

The problem of minimizing the supersonic drag for a given 
lift on a fixed plan form has been approached in different 
ways. Jones, in references 1 and 2, makes use of reverse- 
flow theorems to derive several simple properties of the 
optimum load distribution and to present as well the optimum 
distribution for elliptic plan forms. Graham, in reference 3, 
shows how the proper use of orthogonal loadings can reduce 
the drag at fixed lift. Orthogonal loadings are loadings of 
zero interference drag. The interference drag between two 
loadings is the total drag of each in the down wash field of the 
other. In reference 4, theorems concerning orthogonality 
and reverse flow are developed, whereas in references 5 and 6 
numerical examples of drag reduction by use of orthogonal 
loadings are given. For delta wings with conical camber the 
optimum shapes are derived by Ritz’s method in reference 7. 

In this report Lagrange’s method of undetermined multi- 
pliers is applied to the problem of properly combining load- 
ings for the least drag at a given lift. By use of this method 
a simply expressed property of the optimum loading is found 
which is an integral form of a property established by Jones 
in reference 1 for reversible flows. Jones’ property of the 
optimum loading is that the downwash on the plan form is 
constant in the combined forward- and reverse-flow fields. 
The best combination of four types of nonsingular loading on 
a delta wing is calculated as a numerical example of the use 
of the method. 


SYMBOLS 

A loading strength parameter 

b span 

c local chord 

C D drag coefficient 

C Dli drag coefficient of ith loading 

Cd,h drag coefficient of interference between ?'th and jth 
component loadings 
C L lift coefficient 

C L ,t lift coefficient of ith loading 

C p lifting pressure coefficient 

M Mach number 

m tangent of semiapex angle 

N number of loadings 

n=/3m sweepback-speed parameter 

R functions of d and n (see appendix) 

S' wing area 

X', Y' arbitrary Cartesian coordinates 
l loading on an arbitrary line, jc n dX' 

s, t integers 

x, y, 2 Cartesian coordinates of lifting surface (see fig. 2) 
a local angle of attack of lifting surface 

0=VM 2 - 1 

e small positive number 

e=- y - 

mx 

X Lagrange’s multiplier 

r plan form 

1 — y root chord of arrow wing 

Subscripts: 

i, j tth, jth loading component 

M minimum among all loadings 

0 minimum among N loadings 

X arbitrary loading 

ANALYSIS 

THEORY 

Consider a superposition of N loadings of the form 

C p =A I C p ,iAA 2 C p ,2-f-A 3 C 1 , ) 3-j- . . .-\-A N C pt .\! (1) 


Supersedes NACA Technical Note 3533 by Frederick C. Grant, 1055. 
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where A is the strength parameter and C p is the resultant 
lifting pressure coefficient at a point on the plan form. The 
corresponding local angle of attack may be written as 

a =^4 I ai+jd 2 a2 + ^i3a3 + . • .-\-A N a N (2) 

The local drag coefficient C P a is a quadratic in A which 
may be integrated over the plan form r to give the drag 
coefficient of the wing. Thrust-loaded singularities at the 
leading edge are therefore excluded from the drag. This 
exclusion is merely for convenience and is not necessary. A 
formula for the drag coefficient is 

6'„=4 f C„a dS=\ S £ C D ,ijA t Aj (3) 

O J T ^ i— 1 j=l 

where 

C'd, ij~~ C D , = (Cp, i a j+ Cv, j a i)dS 

The average lifting pressure coefficient on the plan form 
is the lift coefficient, which is 

Cl—-q r C„ dS= y Cl, {A i (4) 

O Jr 1=1 

The problem is to find the set of A’s which yields the 
minimum value of C D subject to the condition that C L is 
constant. Because of the quadratic nature of C D and the 
linear form of C L , Lagrange’s method of undetermined mul- 
tipliers is particularly suitable for the solution as it leads to 
a set of linear algebraic equations. 

As shown in reference 8, a function of the A coefficients 
F = C D +^C L is formed, where X is Lagrange’s multiplier. 
The minimum value of F as determined by the N linear 
dF 

algebraic equations =0 plus condition (4) is Lagrange’s 

solution. In matrix form these equations are: 


2 Cd, i 

Cd, 12 

Cd, 13 

Cd, in 

Cl, i 


~Lir 


~0~ 

Cd, 12 

2 Cd, 2 

Cd, 23 

C D, 2 N 

C L , 2 


-42 


0 

C D , i3 

C-D. 23 

2 C’d, 3 

Cd, 3 n 

Cl, 3 


Az 

_ 

0 

Cd, in 

Cd, 2n 

Co, 3 N ■ 

‘ZCd.N 

Cl,N 


An 


0 

-Cl, i 

Cl. 2 

CO 

• Cl.N 

0 _ 


_ X _ 


Cl. 


(5) 


The equations may be written more simply if first the 
interference drag between the optimum loading and the fth 
component of the loading is computed. From equations (1) 
and (2), the following expressions may be written: 


C Pj o“i— AiCp, iOii J rA2C p _2 OL i J rA i C Pt iCii J r ...+ 

A i C p , i a i “f" . . . A.nC p ,N®-i 

a oC p , t=Aia 1 C v , i-f-AotyCp, i~hA 3 a 3 C Pt ,+ . . . + 

Ajpi^C p ,i~ j~ . . • -j- An&nC p , i ^ 


(6) 


Adding equations (6) and integrating over the plan forrr 
gives 


C D ,0i~-g J" (C p ,«a i +a 0 C p , t )dS=A l C D , uAA.Cn, 2i + 

A 3 C D , 3i+ . . . ~\-2AiC D , i+ . . . -\-A N C D ,m 

= 2D AjC d , n (7) 

1=1 

This expression for C Dl0i is a part of the left-hand side of 
the zth equation of the linear set which is now written as 

C D ,oi+^C L ,i=0 (8) 

A simple property of the optimum load distribution may 
now be derived. First C D , o is rewritten by use of equation (7) : 

Co,o—n A ( Cd, oe (9) 

& !=1 

or using equations (8) and (4) 

C D .,= -\\C L (10) 

Substituting equation (10) into equation (8) gives 

Cd, oi = 2 ,/ Cl, i (11) 

t-r, 

Since equation (11) holds for any number of loadings, let 
the number of components increase without limit to include 
all possible loadings. For an arbitrary loading X and the 
absolute minimum M, equation (11) may be written as 

n ..... n Cd, m 
Cd,: mx — = ~7j '-l.x 

Ul 


( 12 ) 
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The meaning of equation (12) may be simply expressed as 
follows : The interference drag between the optimum loading 
and any loading at the same lift coefficient is constant. If 
the reversibility theorem is applicable, equation (12) is an 
integral equivalent of a condition established by Jones in 
reference 1. Jones’ condition states that for the optimum 
loading the downwash on the plan form is constant in the 
combined forward- and reverse-flow fields. Barred variables 
will represent the reverse flow which has the same lift load- 
ing on the plan form but, in general, a different surface 
shape. Then, by reversibility, 

J 0 T , M a x dS= J <%, ,,a Y d-S=J’ a,,C v , x dS (13) 
By definition, C D<MX is 

Cl)'MX — -g ( ( 'p, M a x “b “v x)dS 

Therefore, equation (12) may be written as 

j t,v.Y(«,»+s»y*s=2 r, ( Y' j (\ s ds (i4) 

Since C,,,x is arbitrary, a M -j- a M must be constant. Hence, 
« U +*M = 2 (15) 

'-i 

This is the condition derived by Jones in reference 1 . Equa- 
tion (12) is then an equivalent integral form of this condition. 

Equation (12) shows the orthogonality of the optimum 
loading to, and only to, zero lift loadings. This result, 
which was stated by Graham in reference 3, is seen to be a 
special case of a more general interference drag property 
given by equation (12). 

COMPARISON WITH THE METHOD OF ORTHOGONAL LOADINGS 

If two loadings are to be combined, it may be shown that 
Graham’s method of orthogonal loadings (ref. 3) and the 
present method are equivalent. If the resultant combina- 
tion of two loadings is combined by the method of reference 3 
with a third loading, the lift ratio of the first two loadings is 
unchanged in the best combination of the three. If w>2 
loadings are successively combined in the manner of refer- 
ence 3, the first n— 1 loadings are not allowed to adjust their 
relative strengths upon addition of the nth. In the present 
Lagrangian method every loading has equal freedom to 
adjust. For this reason, the Lagrangian method should 
be more rapidly convergent. 


NUMERICAL EXAMPLE 

Tucker in reference 9 presents formulas for the surface 
coordinates of delta and arrow wings which support four 
types of pressure distribution. The formulas are given for 
subsonic leading edges and supersonic trailing edges. In 
the notation of this report (fig. 1) a combination of the four 
loadings may be written : 

C^A.+A.x+A, ^+A t £ (16) 

Formulas for the C D , tj quantities may be derived from equa- 
tion (16) and the surface formulas given in reference 9, by 
integrations over the plan form. Details are given in the 
appendix. 

The optimum-drag results are presented in figure 2 along 
with the corresponding drag values for a flat delta wing 
with and without leading-edge thrust (ref. 10). The drag 
values for the four component loadings taken alone are also 
shown. In addition, the drag of the conically cambered 
optimum delta wing (ref. 7) and Jones’ absolute minimum 
for narrow wings (ref. 1) are plotted. The optimum A 
values are tabulated in the appendix. 

Noteworthy in figure 2 is the closeness with which all the 
optimum drags agree with each other and with the drag of a 
flat delta wing which has a thrust-loaded leading edge. The 
close approach of the present optimum of four loadings to 
Jones’ absolute minimum for narrow wings is also evident. 
The data indicate that the relatively low drag of the flat 


V 



Figure 1. — Arrow plan form. 
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delta wing with leading-edge thrust can be equalled by prop- 
erly combining a few loadings having finite pressures every- 
where on the plan form. A plausible speculation suggested 
by the data is that it is possible to come very close to the 
minimum drag on a delta wing with but a few steps in a series 
approximation. Perhaps, too, a restricted minimum, such as 
the one for conical camber, gives a close approximation to 
the absolute minimum drag if the restriction is not too 
unnatural. 

Since the vortex drag of a wing at any Mach number 
depends only on the spanwise loading, a departure from the 
elliptic spanwise loading is a measure of the vortex drag in 
excess of the least possible drag. In figure 3 the spanwise 
loading of the optimum combination is shown at the ex- 
tremes of the sweepback range. There is good agreement 
with the elliptic loading especially for the case of extreme 
sweepback (n = 0). Because for extreme sweepback the 
wave drag vanishes, a direct comparison of the vortex drag 
of the optimum combination and the elliptic spanwise load- 
ing is given by figure 2 at w = 0. The elliptic spanwise load- 
ing has the drag parameter value ~ 

It is shown in reference 2 that the wave drag due to lift 

depends on all the lift loadings l(Y') where 1= J* C v dX' and 

X' is an arbitrary direction inclined to the free stream at 
more than the Mach angle. The coordinate Y' is perpen- 
dicular to X'. A sufficient condition for minimum wave 
drag is shown to be that l(Y') is an ellipse. In figure 4 the 
loading of lines perpendicular to the free stream, or chord- 
wise loading, is shown for the optimum combination with a 
sonic leading edge (» = 1 ) . Agreement with the elliptical 
loading is poor. For the case of extreme sweepback (n= 0) 
no chordwise loading for the optimum combination is shown 
in figure 4 because it is partially arbitrary. (See appendix.) 
The allowable variations of the optimum loading at n— 0 
correspond to changes in the oblique loadings that do not 
change the spanwise loading. This result emphasizes the 
vanishing of the wave drag with extreme sweepback. 



Fiocre 3. — The loading of lines parallel to the free stream for the 
optimum combination. m — Ci.— 1. 
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Figure 4. — The loading of lines, perpendicular to the free stream for 
the optimum combination, m = Ci, = 1 . 

CONCLUDING REMARKS 

Lagrange’s method of undetermined multipliers is applied 
to the problem of properly combining lift loadings for the 
least drag at a given lift on supersonic wings. 


The method shows the interference drag between the 
optimum loading and any loading at the same lift coefficient 
to be constant. This is an integral form of the criterion 
established by Robert T. Jones for optimum loadings. 

The best combination of four loadings on a delta wing 
with subsonic leading edges is calculated as a numerical 
example. The loadings considered have finite pressures 
everywhere on the plan form. At each Mach number the 
optimum combination of these four nonsingular loadings has 
nearly the same drag coefficient as a flat plate with leading- 
edge thrust. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va„ July 27, 1955. 




APPENDIX 


DETAILS OF NUMERICAL EXAMPLE 


INTERFERENCE DRAG FORMULAS 


Inasmuch as the pressure coefficient C p ancl corresponding 
angle of attack a are given by 


\V\-LA 


C„=A,+A,x^A, \ m +A 4 — 2 I 

a— A\Ot\-\- AvOLi-^r AiOti J 


(At) 


then the local drag coefficient mat' be written as 

C p cx=Ai 2 (oi 1 )-j r AiA2(xoc l -\-a 2 )A-A ] A 3 «i+« 3 ^+ 


)+ J / l 2 2 (r« 2 )+yl 2 /1 3 

(i+#) 2 +2(i — ft 2 )(0+e 2 ) 

1 +ft 2 0 

/ \m / 

2 a 1— ft 2 

k’lTffj 


A- A 


{i a * +xa ) +A ^t a ) + 

AaAi (h “ 3+ m (jn 2 “ 4 ) 


(A2) 


The, required C D , tj functions are the averages over the plan 
form (fig. 1) of the quantities in parentheses in equation (A2). 
Rather than a, itself, reference 9 gives the surface ordinate 
2 ,- which is the chordwise integrated value of a,-: 


h=— |*«f 


dx 


(Ad) 


The values given for r,- are 


- /> 
m 1 


„=— II, 
m 


z 3 =— Il s 
3 m 


‘ Ci m 4 


(A4) 


V 

The values of R t are functions of d = — tabulated in refer 


mx 


■once 


9 for different values of n. The equations for R, arc 


R\=^ jjih 1— ft 2 # 2 — 2 cosh 1 
a^I— ft 2 (l+#) cosh -1 


1 +n 2 e 
n(l+d ) 


+ 

+ 


Vl— ft 2 (l — 6) cosh 1 


l-ft 2 # ~ 

a(I —e) _ 


(A5a) 


7?,= — 1 —ri 2 d' 2 —2d 2 cosh 
1 [~n 2 (l-0 2 ) 

A 1— n 2 L 2 
1_ rn 2 (l-e 2 ) 

a'1 —nl L 2 


+ 


-T#-f,# 2 J cosh 1 
— #+0 2 ~] cosh -1 


l+ft 2 0 

n{\+6) 


+ 


_1 ~n 2 e I \ 

ft(l-#)l J 


(A5b) 


R 3 = -, 1 fj- Vl 1 + 3# 2 -J ft 2 # 2 ) cosh - 1 


47 r \_2 


+ 


{\~ey-2{\~n 2 )(e-e 2 ) 


cosh -1 


i-rfer\ 

<i-#)U 


n* 


f(i-ft 2 e 2 )^ 

1 ° 1 fYn 2 

li, JU/I, 2o2/-i „2a2. aoi 

i 

\ 3(1 -ft 2 ) 

| ^ o/h i-\\ t L U \ X l v U Uv vUoll 

3n(l— tt) : 

ft# 


(A5c) 


+ 


(1 


1 P 

— ft 2 )- L 


'6— 9ft 2 +2n 4 


9 Q r) 2 

(# 2 +# 3 ) + g ^ (#-# 3 )- 


6 

2- 3a 2 


(1+# 3 ) cosh 1 


\_ + n 2 e 

«(!+#) 


+ 


(1 


i n 

-n 2 )’» _ 


6— 9ft 2 +2ft 4 


0 # 2 -# 3 )- 


2 ' 3 " (d-0' s )—j (1~# 3 )J cosh 1 (A5d) 

For terms in equation (A2) of the type (? ijm) s a u a spamvise 
integration of z, gives the following average on the plan 
form : 

U ,Vi[Va- 


( 1 - '‘ ) ' +,+, raw - rf “] ,A6) 


For terms of the type xa. t an additional integration by 
parts in the x direction is required to maintain the R t func- 
tions intact under the integral signs. The result for this 
case is 


i f xotf ds= 2 - - r^ ) -a-M) ,+2 f t,— 

S Jr l- f l«L< + 2 Jo (1 — (0 


0 # ) 


y #)'+ 3 dd+ 


(i- l i) ,+, r /?*(#) 


£+2 


< A7) 
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In formulas (A6) and (A7) the value of t for each i Is as 
follows : 


i 

t 

1 

1 

2 

2 

3 

2 

4 

3 


By applying formulas (A6) and (A7) to the integration of 
(A2), the following equations for C D>ti are derived: 

2mC DA =~^ i? l( l)-4(l-,)j;^ dd (A8a) 

2(i -'‘ ), X'<uSi M <Asb) 


(i-^) 

6L\(e) 






2 <‘-*>x:<£ 


(?L 

M®) 


d6» 


toCd 




2 mC D , 


1 2(1 
2 ( 1-m ) 3 J o ‘ 
*»(1) 


M0) 5 


(1-K0) 1 




(A8c) 


tf0- 


(A8d) 


1 M 


4(1 -m) 3 


mCj 


J„(1- M 0) 6 ^ +(1 M) Jo(l-M0) 4d0 

(A8e) 

■"■“=2(Uri [ft(i)+A<i)]-2(i-,)«/; 4 .- 

2<i -' ), X'(^* + 5 (i -^X , (w^ <as,) 


2 < i -“) , X(uS? *+* a ->x: cuS?'* 


(A8g) 


(ASM 


1 0 2 ff 3 (fl) 

M*) 6 


r/0- 


mCo ' 34= 5(l^ ) [«3(l)+i?4(l)]-2(l-M) 4 J o 1 ~ 

w-rf! 0%* < Asi > 

2m<7 »'*-3 (r w-Aa-crX’ PS * < A8 >> 

The required C Lii functions are simple integrals over the 
plan form which yield 


C LA =\ 


2— m 


C,. 2 = 


t 7 A 

L - 3 3 


n i 

° i,4_ 6 


NUMERICAL CALCULATIONS 


(A9) 


The integrals in equations (A8) were, in general, evaluated 
numerically. However, several of the integrands in equa- 
tions (A8) have the form and These 

(1— ju0)‘ (1 

functions have an infinite discontinuity at 0=0. For such 
a discontinuity, numerical methods break down. Near 
zero the following approximation is integrated analytically : 


K l (0)«if 1 (e)+^ cosh- ~^cosh-» i 
cosh- ^+^cosh- 1 


Xo<e^€«:i (Aio) 


The integrals for the region O^0^e can then he ap 
proximated : 


•^ > 3 I 

Jo (1- 


M 

2 TV 




(All) 


where 

/(,)- X' 


d0 


(1- M 0)‘ 


H 


e , t(t + l)in 2 , 
^ 2^ 2! 

<(f+l)(f + 2) Jl 
3! 4 


+ 


-] 


and 


cosh 1 -i; 
n8 


J t LUOU 

.7I=S 


d0 


(A 12) 


(A 13) 
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The integral in equation (A13) may be evaluated by 
expanding the denominator by the binomial theorem and 
writing /(e) as an infinite series 


TTiOd. I3 : 


4 4 1 , 

3 3 loge 
4ir 


(A 18c) 


where 


= aolO + Ulfl + ®2*2 _ b®3f3 _ l _ • • • 

(A14) 

Oo=T 

(A 15a) 

t 

«1=Y M 

(A 15b) 

t(f + l) ■> 
<h= 2T M- 

(A 1 5c) 

f(t + l)(f + 2) 3 

a 3 g, M 

(Al5d) 

t(t + l) . . . (t+s— 1) 
s s\ M 

(A 1 5e) 

i 0 = 1 cosh -1 ~ dfl 
Jo nd 

(Al5f) 

ii— \ d cosh -1 — - dd 

Jo nd 

(A 1 5g) 

u= ( d 2 cosh -1 \ dd 
' Jo nd 

(Al5h) 

i s = 1 6 s cosli -1 4 dg 
Jo nd 

(A 1 5i) 


The i s integrals of equations (A15) are evaluated by use 
of the relation 

r i fl s+1 i i r f) s 

0 s cosh- 1 -XT cosh ~ «+oJLi l=¥% de ( Al6 ) 

J nd s+1 nO s+1 J H— n 2 d 2 


EXACT CALCULATIONS 


mC D , i 


4, . „ 5 
3 l0ge 2 6 

4ir 


2m,C D ,2= -j- 

47 r 


mC L 


4 4 , » 

3“3 l0S ' 2 


4 x 




2 mC DA 


4 4, 0 

3~3 l0ge2 
4ir 


mCj 


D, 34- 


80+15 1, >gc 2 

47T 


2 mCi 


D, 4- 


4ir 


(Al8d) 

(Al8e) 

(Al8f) 

(A 1 8g) 

(Al8h) 

(Al8i) 

(AlSj) 


Equations (18) provide the interference drag coefficients re- 
quired to calculate the vortex drag due to any combination 
of A values. 

In the solution for the optimum A values, the parameters 
Ai, A a , and A :s are found to be linearly related and one of 
them may therefore be chosen arbitrarily. Choosing A\ 
yields : 


At the extremes of the sweepback range, equations (A8) 


may be evaluated exactly, 
back (n=0), there results: 

For the 

case of extreme 

sweep- 

R '=-L( 2+los ° i^e 2+e loge l+£ 

) 

(Al7a) 


(l+0 loge \-^r e + e ' 2 lo ?e 1 

- e v 

(A 17 b) 


[| + (l+30 2 ) loge 

e 4 

jl+ 

CL 

c 

(A 17c) 

■Vi-0 2 ^ 


;+40 2 +(0+20 3 ) loge ^4 

-3d 2 loge , 

(Al7d) 


2mC Dt i= 

loge 2 

7 r 


(A 18a) 


771 Od, 12” 

2.4, 
3+3 l0g 

e 2 

(A 18b) 



4r 



Ai=A { 

(A 19a) 

. 3(40-1) 

1 +a~ 2Ai 

(A 19b) 

, , 3(38a— 26a 2 — 11) 

31 (2-3a)(l+o) 

(A 19c) 

30 (3a — a' 2 — 1) 
J 4_ (2-3a)(l+a) 

(Al9d) 

9 (4a — 1 ) (3 — 2a) ( 1 — 2a) 
D -°~ 8tt (2— 3a)(l+a) 

where 

(Al9e) 

a=| (1 loge 2) 


The spanwise loading may be written as 


^-rr+A^+2/ 2 (A 4 )J (1 — y) 

?/> 0 (A20) 




when m=C L = 1 . Substitution of the A values given in 
equations (A19) shows that l(y) is independent of the vari- 
ations in Ai > ^- 2 , and A 3 . 

For the case of a sonic leading edge (n=l), 

( 271^—2 cosh” 1 (A21a) 

R 2 =—~(^2-yJl-e 2 3 4 5 —2d 2 cosh- 1 (A21b) 

^—cG^'-O+H coslr ‘ 0 <A21c) 

"-“sK+f” 1 ) “ slr ‘ ?] < A21d > 

2mC D ,i =\ (A22a) 

mC D , 12 =| (A22b) 

mC D , 13=^+^ (A22 c) 

mC D . u =jQ (A22d) 

2 mC D ,t=\ (A22e) 



(A22f) 

mC * M== m 

(A22g) 

2mC D . ? =~ 

(A22h) 

mC D , 3i =±+ JL 

(A22i) 

2mC D , 4 =ggQ 

(A22j) 


CALCULATED VALUES OF A 

The table that follows contains the calculated values of A 
for the optimum combination through the sweepback range. 
Four significant figures are given, since the tabidated values 
of R have four decimals. Values of C D , 0 for m=C L = 1 are 
also shown : 


— 

, 



- - — 

— 


-i, 

a 2 

*' l 3 

•i. 

c 

K 'U,0 

0 

.2 

1.993 

-2. 580 

1.435 

1. 054 
1.517 

0. 0830 
. 0899 

.4 

1.977 

-2. 571 

1.590 

1.244 

. 1105 

.6 

1.781 

-2. 147 

1.472 

.9568 

. 1398 

.8 

1.G41 

-1.818 

1.304 

.6957 

. 1700 

1.0 

1.357 

-1.201 

1. 259 

. 1406 

. 2295 
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